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ABSTRACT 

We present the first mid-infrared spectra of brown dwarfs, together with observations of a low-mass star. 
Our targets are the M3.5 dwarf GJ 1001 A, the L8 dwarf DENIS-P J0255-4700, and the T1/T6 binary system e 
Indi Ba/Bb. As expected, the mid-infrared spectral morphology of these objects changes rapidly with spectral 
class due to the changes in atmospheric chemistry resulting from their differing effective temperatures and 
atmospheric structures. By taking advantage of the unprecedented sensitivity of the Infrared Spectrograph on 
the Spitzer Space Telescope we have detected the 7.8 /im methane and 10 /mi ammonia bands for the first time 
in brown dwarf spectra. 
Subject headings: stars: low mass, brown dwarfs - infrared: stars 



1. INTRODUCTION 

Within the past few years the number of known brown 
dwarfs has increased dramatically, with both the Sloan Dig- 
ital Sky Survey (SDSS; Yorketal. 2000) and Two Micron 
All-Sky Survey (2MASS; Sfa-uts kieetalJll997l) accounting 
for the bulk of the discoveries. These objects have been iden- 
tified based on broad-band photometry with follow-up opti- 
cal and near-infrared spectroscopy. The integrated data from 
these investigations has yielded important information about 
the variation of the chemical composition and of the effects of 
condensates in brown dwarfs in term of their spectral types (e. 
g. Kirkpatrick et al. 2000; Leaaett et al. 2001; Burgasser et al. 
2002t lMcLean et all2003 : Knap p et all2004l) . 

Spectroscopic observations in the mi d-infrared yield im- 
portan t additional information (Saumon, Ma rlev. & Loddersl 
2003), but due to the very low infrared brightness of these 
objects and the effects of the Earth's atmospheric absorption 
these measurements have been difficult or impossible from 
ground-based observatories. With the launch of the Spitzer 
Space Telescope (Werner et al. 2004) and the onboard In- 
frared Spectrograph (IRS; Houck et al. 2004) on 2003 August 
25, we now have an observatory with sufficient sensitivity to 
undertake these studies. Accordingly, the Science Team of the 
IRS has organized a comprehensive observational program of 
M, L, and T dwarfs, covering the range from the earliest M 
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dwarfs to the latest T dwarfs. We present here the first results 
from this program: observations of the M3.5 star GJ 1001A 
(LHS 102A), the L8 dwarf DENIS-P J0255-4700 (hereafter 
J0255-4700), and the T1/T6 binary system e Indi Ba/Bb. 

2. OBSERVATIONS AND FLIGHT DATA REDUCTION 

The three objects described here were observed with the 
IRS instrument on the Spitzer Space Telescope as part of one 
of the major IRS Science Team Guaranteed Time Observing 
programs. The IRS is capable of low-and moderate-resolution 
spectroscopy, with only data from the low-resolution modules 
being reported here. Each of the IRS low -resolution modules 
is fed by two long slits, one used for first order and the second 
used for second order. The IRS low-resolution modules have 
slit lengths of 54. "6 and 151. "3 for both orders in the short-low 
(SL) and long-low (LL) modules, respectively. The spectral 
resolution of the SL module as determined by the slit width is 
AA = 0.06 /im for second order and AA = 0.12 /im for first 
order. The spectral resolution of the LL module as determined 
by the slit width is AA = 0.17 /mi for second order and AA = 
0.32 /im for first order. The slit widths of both of the IRS low- 
resolution modules a re Nvquist sa mpled by the pixel pitch of 
their infrared arrays. iHouck et alJ (J2004) describe the design, 
performance, and operation of the IRS in more detail. 

Table 1 provides a summary of the observations reported 
here. As all of these objects have relatively large proper mo- 
tions, the IRS blue peakup array was used to locate the current 
positions of these objects and to reposition them into the IRS 
slits with an estimated radial uncertainty of approximately 
0"4, (la). Each object was observed at two locations in each 
of the IRS low-resolution slits, with the positions located at 
approximately the 1/3 and 2/3 positions along the long axis of 
the slits. 

The standard IRS data reduction pipeline version S9.1 at 
the Spitzer Science Center processed the data. Stellar stan- 
dards were also observed to calibrate the spectra and re- 
move any spectral artifacts which may have survived the pro- 
cessing; the results reported here were corrected with the A 
dwarf standard star <5 UMi. The zodiacal and galactic back- 
grounds were subtracted before the spectral extractions using 
observations when the target was located in the alternate slit. 
Finally, the spectra were extracted from the sky-subtracted 
two-dimensional array images using an early version of the 
SMART IRS data-reduction software package being devel- 



TABLE 1 

Table of Observations 



Object 



Observation Dates 
Integration Time 



Spitzer AOR Key 



GJ 1001A 



DENIS-P J0255-4700 



e Indi Ba/Bb 



2003 December 16 

SL2 480 seconds 
SL1 480 seconds 
LL2 480 seconds 
LL1 960 seconds 

2004 January 7 

SL2 480 seconds 
SL1 480 seconds 
LL2 480 seconds 
LL1 960 seconds 

2003 September 26 

SL2 56 seconds 
SL1 56 seconds 



4190464 



4192000 



6625792* 



"Observations of the b inary system e Indi Ba/Bb i Scholz et al. 2003; 
IMcCaughrean et all2 004) were taken during the Spitzer In-Orbit Checkout 
and Science Validation period which accounts for it being observed only by 
the IRS Short-Low module and for such a short integration time. Additional 
observations of this object will take place later in the Spitzer mission. 



oped at Cornell University (Higdon et al. 2004). 

Figure 1 shows the resulting IRS low-resolution spectra 
for the three objects reported here. At this early stage in 
the Spitzer mission, our estimated uncertainty in the abso- 
lute photometry is ±25%, but the uncertainty in the spectral 
shape is much less than this. The caption for Figure 1 gives 
the uncertainties in the spectral shapes which were estimated 
from a combination of the measured noise in the individual 
reads of the pixels in the detector arrays and the variation be- 
tween the spectra extracted from the observations at the two 
slit positions. These uncertainties therefore incorporate some 
measure of the residual systematic uncertainties and spectral 
flat-fielding errors. Uncertainties in the absolute flux calibra- 
tion were estimated in different ways for the three objects and 
are described below. 

3. DISCUSSION 

As the atmospheres of brown dwarfs cool with time, their 
spectral signatures reflect a progression of changes in their 
atmospheric chemical equilibrium and condensate structure. 
The three dwarfs discussed here provide snapshots of this 
progression. In an M dwarf such as GJ 1001 A, the ele- 
ments O, C, and N are predominantly found in H2O, CO, 
and N2 a nd the atmosphere is too warm for c ondensation 
of solids lAllard & Hauschildll 119951) : lLoddersI 119991) 1. As 
the effective temperature (T e ff) falls, a variety of conden- 
sates form in the atmosphere, most notably iron and sil- 
icates. These condensates are not expected to be well- 
mixed through the atmosphere, but rather be found in rel- 
atively thin, discrete cloud layers overlying the condensa- 
tion level (Ackerman & Marlev (2001); Marlev et al. (2002); 
iTsuiil ( 120021) : IWoitke & H elling (2004)). As the T eff falls to 
that of a late L dwarf, such as J0255-4700, the cloud layer 
is optically thick and affects either directly (as a major opac- 
ity source) or indirectly (by altering the atmospheric temper- 
ature/pressure profile) all spectral regions. The exact spectral 
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FIG. 1 .— The observed IRS spectra of GJ 1001A, J0255-4700, and the e 
Indi Ba/Bb system. The fluxes for e Indi Ba/Bb have been divided by a factor 
of 10 for clarity. The average la errors in the data points are as follows: for 
GJ 1001A, 2.60 mJy from 5.3 to 14.1 /im and 1.61 mjy from 14.1 to 35.0 
/im; for J0255-4700 they are 0.73 mjy from 5.3 to 14.1 fim and 1.00 mJy 
from 14.1 to 35.0 (im; for the e Indi Ba/Bb system they are 2.26 mjy from 
5.3 to 15.3 fim. 



signature of the cloud depends both on its vertical thickness 
and the particle size distribution of the condensates. In ad- 
dition as the atmosphere cools chemical equilibrium begins 
to fav or first C H4 over C O and then NH3 over N2 (Tsuii 
( 1964); Feglev & Lodders ( 1996); Burrows & Share ( 1999); 
Lodders (1999); Lodders & Feglev (2002): iBurrows et alJ 
(2001)). The strong mid-infrared bands of these molecules 
begin to challenge the domination of water opacity, at least 
in some spectral regions, by the late L's. By the early to 
mid T dwarfs, such as the e Indi pair, the condensate cloud is 
forming quite deep in the atmosphere. In the relatively clear, 
cool atmosphere above the cloud, chemical equilibrium be- 
gins to strongly favor CH4 and NH3 and their spectral fea- 
tures, along with particularly strong bands of water , grow 
in prominence (Marlev et al. (1996); Burrows et al. (1997); 
Allard et ail (120011) : IBurrows. S udarskv. & Lunind 120031) 1. 

Figure 1 shows the resulting pronounced changes in the 
mid-infrared spectral morphology between the three objects. 
The spectra and comparisons with synthetic spectra are con- 
sidered for each of the objects in turn. 

3.1. GJ1001A 

GJ 1001 A is classified as an M3.5 V star 
(Hawl evTCiizis. & Reid I 119961) and is part of a binar y 
system with the L5 dwarf GJ 1001B (lGoldmanetal.ll 19991) . 
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FIG. 2.— The observed IRS spectra of GJ 1001A, J0255-4700, and the e 
Indi Ba/Bb system compared with model spectra. The fluxes for e Indi Ba/Bb 
have been divided by a factor of 10 for clarity. The observed data are denoted 
by the solid black lines while the model spectra are in the thicker blue lines. 
The various model parameters for each of the spectra are given in the text. 
The most noticeable HiO, CH4, and NH3 bands are marked, although weaker 
H2O absorption extends from the blue edge of the figure to 7.0 fim, and NH3 
absorption in e Indi Ba/Bb extends from 8.5 fj,m to the red edge of the figure. 
In the case of GJ 1001A the entire displayed spectrum is dominated by water 
absorption, while in J0255— 4700 water is dominant longward of 9 11m. In 
a Tl (the brightest component of e Indi B), NH3 is present but not strong 
enough to be visible in the 12-15 /im region, so that longward of 12 /im the 
spectral morphology is still primarily due to water absorption. The strong 7.8 
fim absorption band from CH4 is evident in both of the cooler objects. 



The distance between the primary and secondary is large 
enough (18") that there is no contamination of the primary's 
spectrum by stray light from the secondary. Using optical and 
near-infrared spectra, Leggett et al. (2002) find a best fit to 
this object with T e ff = 3200 K and log g = 5 using a mixing 
length parameter of 2 in the AMES-dusty models of Allard 
et al. (2001). Figure 2 compares this same model, smoothed 
to a spectral resolution of R ~ 400, with the IRS data. The 
6.5 /im feature seen in the IRS spectrum also appears in the 
model and arises from H2O opacity. 

For the o bject GJ 1001A the near-infrared observations 
(Leaaett et al. 2002) can be used with the model and the re- 
sulting predicted photometric flux compared with our obser- 
vations. Although the shape of the model spectrum is an ex- 
cellent match to the IRS data, the model fluxes derived from 
the near-infrared observations had to be scaled by a factor of 
0.87 to fit the observed data. The near-infrared flux levels are 
consistent with the T e jj and radius derived from the luminos- 
ity if the dwarf is older than 0.15 Gyr (Leaaett et al. 2002; 
lBaraffeetalJl9 98). 



3.2. J0255-4700 

J0255-4700is classified optically as an L8 dwarf. 13 Unfor- 
tunately the trigonometric parallax distance of this object has 
not been measured. Considerable structure can be seen in the 
IRS spectrum and to interpret these resul ts we have compared 
the IRS data with the model described bv lMarlev et all (120021) 
and Saumon et al. (2003), with T eff = 1400 K, log g = 4.5 
in Figure 2. In this model atmospheric condensate clouds are 
treated as described by Ackerman & Marlev (2001) with a 
sedimentation efficiency parameter f se( j = 3. The model has 
been normalized to match the observed flux in the 8.6 - 9.0 
^m region. Models that include condensate sedimentation 
in the chemical equilibrium but ignore cloud opacity require 
unrealistically high T e ff to acceptably fit the data and these 
are not shown here. We judge the excellent agreement of 
the cloudy model to the J0255-4700 data as further evidence 
for finite-thickness silicate and iron clouds in the observable 
atmospheres of late L dwarfs. As can be seen in the figure 
the correlation between the model spectrum and the IRS data 
is excellent, although there is a slight discrepancy between 
the data and model predictions in the H2O feature at 5.5 - 
6.5 /im. In addition there is a broad shallow dip in the IRS 
data compared with the model between 9 and 11 /im that 
may be due to a silicate cloud feature. The IRS spectrum 
shows the first detection of the 7.8 /im CH 4 band in a brown 
dwarf. CH 4 was detected at 3.3 /im in L dwarfs by No ll et alJ 
(2000) and since the 7.8 /im band is of comparable strength 
to the 3.3 /im fundamental, it is not surprising to find this 
mid-infrared signature of atmospheric methane in an L dwarf. 
Mo dels with (Saumon, Marlev, & Lodders (2003)) and with- 
out dBurrows et alJ (11997b : lAllard et aTT (2001 V) dust opacity 
also anticipated this feature in the L dwarf T e ff range along 
with the water features noted above. Methane absorption is 
even more obvious in the spectrum of the e Indi B system dis- 
c ussed below. 

ICreech-Eakman et al.ll|2004l) used OSCIR at Gemini South 
to image J0255-4700 in the JV-band (~ 8.1 - 13.4 /im) and 
three narrowband filters centered at 8.8, 10.3 and 11.7 jum. 
Flux density detections were obtained in the N and 8.8 /mi 
bands, with upper limits for the 10.3 and 1 1.7 /mi bands. De- 
spite the estimated 25% uncertainty in our absolute flux cal- 
ibration, the IRS spectrum cannot be scaled to match both 
their A/-band and their 8.8 /mi flux density measurements. 
Creech-Eakman et al. (2004) note that the N-band flux den- 
sity is a factor of 2 to 3 greater than the narrowband values 
and discuss possible causes of this discrepancy. Our data can 
rule out two of their suggestions - the presence of peculiar 
emission or absorption lines. Figures 1 and 2 show that there 
are no strong emission features in the N bandpass and also 
that the narrow bandpasses do not fall in regions of strong 
absorption. While it is difficult to reconcile all of the Creech- 
Eakman et al. photometry with our IRS spectrum, our data 
can be brought into approximate agreement with the three nar- 
rowband measurements if the spectrum is scaled to the lower 
limit of the 8.8 /im OSCIR detection. This may indicate that 
the A/-band measurement is in error. 

3.3. e Indi Ba and Bb 

The e Indi Ba and Bb system consists of a Tl and T6 
dwarf separated by 0."732 (McCaughreanet al. 2004). The 
IRS does not have the spatial resolution to separate the two 

1 3 http ://spider. ipac . cal tech. edu/staff/davy/ ARCHIVE 



objects so the measured spectrum is a composite of the 
two individual spectra. In this object we find the first evi- 
dence for NH3 absorption in the spectrum between 10 and 
1 1 ^m, although its presence in very cool brow n dwarf at- 
mospheres has been expected for some time (Tsuii 1964; 
Marlev et al. 1996; Fealev & Lodders 1996; Saumon et al. 
2000; Allard et al. 2001; Burrows, Sud arskv. & Lu nine 2003; 
Saumon, Marlev, & Lodders 2003). Saumon et al. (2000) re- 



ported a weak detection of NH3 in the near infrared spectrum 
of the T6 dwarf Gl 229B but this is the first strong detec- 
tion of NH3 in a brown dwarf. With the lower T e yy for the 
components of this binary system the condensate clouds are 
predicted to be below the photosphere and thus have a less 
pronounced effect on the composite spectrum than they do on 
the L8 dwarf's spectrum. 

A composite model spectrum provides an excellent match 
to the IRS data and was obtained as follows. As with 
the model for J0255-4700 we used the model described by 
iMarlev et al . (2002) and Saumon et al. (2003). From the 
estimates of Li, i for each component (McCaughreanet al. 
2004) and the estim ated age of e Indi A (0.8-2 Gyr; 
Lachaume et al. 1999), we use our brown dwarf cooling cal- 
culations (Marlev et al. 2004) to obtain (T eff , log g, R/Rq) = 
(1250 K, 5.13, 0.094), (840 K, 4.89, 0.100) for the Tl and the 
T6 dwarfs, respectively. Selecting the closest cloudless syn- 
thetic spectra in our grid of models (1200 K, 5.0) and (800 K, 
5.0) which are within the uncertainties of those parameters, 
and using the above radii and the trigonometric parallax of e 
Indi A (Perrvmanet al. 1997), the absolute flux at the Earth 
of the T dwarf binary was obtained. Note that in Figure 2, 
the composite model spectrum has not been normalized to the 
IRS data, so that in this case the predicted flux agrees with the 
IRS data to within the data errors. 

4. SUMMARY AND CONCLUSIONS 

(1) We have taken advantage of the unprecedented sensi- 
tivity of the Infrared Spectrograph instrument on the Spitzer 



Space telescope to observe two brown dwarfs. These are the 
first spectra in this wavelength range ever reported for this 
type of object. The spectrum of an M dwarf is also reported 
here, providing a contrast between objects with masses above 
and below the Hydrogen-Burning Minimum Mass. 

(2) The observed mid-infrared spectral morphologies of 
these three objects vary strongly with their spectral classes, 
reflecting the changes in gas chemistry driven by the temper- 
ature of their atmospheres. 

(3) Model comparisons with the observed spectra of all 
three objects show good agreement, with only a few minor 
deviations. Models with cloud opacities do a better job of re- 
producing the data from J0255-4700 than do models without 
clouds. 

(4) We report here the first positive detection of the 7.8 
/immethane band. We also report the first unambiguous de- 
tection of NH3 in a brown dwarf atmosphere. 
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